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Abstract 


This  report  presents  an  evaluation  of  simulated  RADARSAT-2  polarimetry  products.  The 
RSAT2SIMU  software  developed  by  MDA  was  used  to  simulate  RADARSAT-2  products  from 
Environment  Canada  (EC)  CV-580  synthetic  aperture  radar  (SAR)  single-look  complex  (SLC) 
data  by  increasing  the  noise  floor  and  degrading  the  image  resolution.  The  evaluation  includes 
assessment  of  the  reduction  in  the  probability  of  missed  detection  for  polarimetric  relative  to 
single  channel  radar  operation  for  ship  detection,  and  the  potential  benefit  of  polarimetric  target 
decomposition  to  generate  ship  target  classification  features  and  to  segment  the  ship  target  of 
interest  from  the  ocean  background.  A  key  recommendation  is  to  develop  algorithms  that 
combine  all  information  available  from  polarimetric  signature  analysis  methods  to  aid  in  ship 
classification.  The  results  show  that  the  ship  detection  and  target  classification  methods 
developed  using  airborne  polarimetric  SAR  data  can  potentially  be  applied  to  RADARSAT-2 
data. 
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Resume 


Le  present  rapport  e  value  des  produits  de  polarimetrie  RADARS AT-2  simules.  Le  logiciel 
RSAT2SIMU  cree  par  MDA  a  servi  a  simuler  des  produits  de  RADARSAT-2  tires  de  donnees 
singulieres  complexes  (SLC)  du  radar  a  synthese  d’ouverture  (RSO)  CV-580  d’Environnement 
Canada  (EC)  par  l’elevation  du  seuil  de  bruit  et  la  degradation  de  la  resolution  d’image. 

L’  evaluation  vise  a  determiner  la  reduction  de  la  probability  de  detection  manquee  de  la 
polarimetrie  comparativement  a  l’utilisation  d’un  radar  monocanal  pour  detecter  des  navires  ainsi 
que  les  avantages  possibles  de  la  decomposition  des  cibles  polarimetriques  pour  generer  les 
caracteristiques  de  classification  des  cibles  de  navire  et  pour  segmenter  une  cible  de  navire 
d’interet  par  rapport  au  fond  oceanique.  Une  recommandation  cle  consiste  a  etablir  des 
algorithmes  combinant  toute  1’ information  obtenue  selon  des  methodes  d’ analyse  des  signatures 
polarimetriques  afin  de  faciliter  la  classification  des  navires.  Les  resultats  revelent  que  les 
methodes  de  detection  des  navires  et  de  classification  des  cibles,  qui  se  fondent  sur  des  donnees 
polarimetriques  de  RSO  aeroporte,  pourraient  eventuellement  etre  appliquees  aux  donnees  de 


RADARSAT-2. 


DRDC  Ottawa  TM  2007-189 


Executive  summary 


Evaluation  of  Simulated  RADARSAT-2  Polarimetry  Products 

Liu,  Chen;  Vachon,  Paris  W.;  DRDC  Ottawa  TM  2007-189;  Defence  R&D  Canada 
-  Ottawa;  September  2007. 

Introduction 

RADARSAT-2,  which  is  expected  to  be  launched  in  late  2007,  will  carry  a  C-band  synthetic 
aperture  radar  (SAR)  that  will  offer  many  operational  modes  such  as  polarimetric  SAR  (PolSAR). 
A  PolSAR  system  alternately  transmits  horizontally  (H)  and  vertically  (V)  polarized 
electromagnetic  pulses,  and  then  measures  both  the  horizontally  (H)  and  vertically  (V)  polarized 
scattered  fields.  There  are  four  possible  combinations  of  incident  and  scattered  electric  fields, 
HH,  HV,  VH,  and  VV.  RADARSAT-2  will  offer  quadrature  polarization  (quad-pol),  selective 
polarization  (HH  and  HV  or  VH  and  VV),  and  selective  single  polarization  (HH  or  HV  or  VH  or 
VV). 

Previous  studies  have  shown  that  there  is  a  significant  improvement  in  ship  detection  sensitivity 
and  a  reduction  in  false  alarm  rate  by  using  PolSAR  data,  especially  for  small  ships.  The  greatest 
improvement  is  obtained  with  quad  polarimetric  data.  However,  the  swath  width  which  can  be 
covered  is  reduced  due  to  system  bandwidth  constraints.  Good  performance  can  also  be  obtained 
with  dual  co-polarization  with  phase,  which  in  principle  provides  a  wider  swath  width  due  to 
reduced  bandwidth  requirements.  Performance  is  characterized  by  the  probability  of  missed 
detection  as  a  function  of  false  alarm  probability. 

The  present  study  investigates  the  feasibility  of  applying  algorithms  for  improving  ship  detection 
and  classification  by  using  airborne  PolSAR  data  converted  to  simulated  RADARSAT-2 
polarimetry  products.  The  PolSAR  ship  detection  algorithm  used  in  this  study  was  developed  at 
DRDC  Ottawa.  Three  target  classification  methods  were  applied:  Pauli,  Cameron,  and  the 
symmetric  scattering  characterization  method  (SSCM).  The  software  for  the  Pauli  and  the 
Cameron  methods  were  developed  at  DRDC  Ottawa,  while  the  SSCM  tool  was  implemented  in  a 
PolSAR  workstation  that  was  developed  at  the  Canada  Centre  for  Remote  Sensing  (CCRS)  and 
was  made  available  to  DRDC  Ottawa  by  the  Canadian  Space  Agency  (CSA). 

The  RADARSAT-2  images  were  simulated  using  the  RADARSAT-2  simulator  (RSAT2SIMU) 
that  was  developed  by  MDA  and  was  also  made  available  to  DRDC  Ottawa  by  CSA.  The 
RSAT2SIMU  software  simulates  RADARSAT-2  products  from  Environment  Canada  (EC)  CV- 
580  C-band  polarimetric  SAR  products  in  single-look  complex  (SLC)  format  by  appropriately 
increasing  the  noise  floor  and  degrading  the  spatial  resolution.  However,  the  simulator  can  not 
account  for  the  specific  geometry  of  RADARSAT-2,  that  is,  the  altitude,  incidence  angle,  and 
swath  coverage.  The  simulated  data  quality  depends  on  the  input  CV-580  SAR  data  quality. 

The  input  CV-580  data  considered  in  this  work  was  obtained  during  the  Maritime  Sensor 
Integration  Experiment  (MARSIE)  trail.  The  data  set  includes  several  known  ships  including  the 
25  m  long  Dominion  Victory. 
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Significance 


This  study  demonstrates  that  the  algorithms  developed  for  improving  ship  detection  and 
classification  using  airborne  PolSAR  data  will  also  be  applicable  to  RADARS AT-2  polarimetric 
data.  The  results  demonstrate  the  significant  improvement  in  ship  detection  performance  that  is 
expected  from  using  polarimetric  data.  Furthermore,  it  is  shown  that  the  polarimetric  signature  of 
a  target  may  also  improve  ship  classification  performance. 


Results 

For  ship  detection,  the  results  show  the  advantage  of  polarimetric  systems.  For  example,  the  dual 
co-polarization  system  with  phase  provides  better  detection  performance  than  single  channel  HH 
polarization.  In  principle,  a  dual  co-polarization  system  could  provide  a  wider  swath  coverage 
than  a  polarimetric  system. 

For  target  classification,  the  Cameron  and  SSCM  target  decomposition  methods  provide  unique 
information  on  the  distribution  and  types  of  scatterers  that  constitute  the  target  of  interest, 
providing  an  approach  for  target  classification.  These  decomposition  methods  work  best  for  well- 
focused,  large  ships.  The  Pauli  decomposition  method  provides  a  simple  but  efficient  way  to 
segment  the  target  of  interest  from  the  surrounding  ocean  clutter. 


Future  Plans 

The  results  demonstrate  that  target  decomposition  and  polarimetric  signature  analysis  may 
improve  target  classification.  To  achieve  target  classification,  the  following  issues  are 
recommended  for  further  study: 

•  The  relationship  between  the  physical  properties  of  the  target  and  the  elemental  scatterers 
derived  from  the  Cameron  and  the  SSCM  polarimetric  decomposition  methods; 

•  The  feasibility  of  applying  polarimetric  decomposition  methods  to  automatic  target 
recognition  (ATR)  algorithms  by  using  scattering  elements  or  estimated  scatterer  orientation 
angles  as  target  features  in  target  classification  software;  and 

•  The  development  of  algorithms  that  fuse  all  information  available  from  polarimetric 
signature  analysis  methods  to  aid  in  ship  classification. 


IV 


DRDC  Ottawa  TM  2007-189 


Sommaire 


Evaluation  of  Simulated  RADARSAT-2  Polarimetry  Products 

Liu,  Chen;  Vachon,  Paris  W.;  DRDC  Ottawa  TM  2007-189;  R  &  D  pour  la  defense 
Canada  -  Ottawa;  September  2007. 

Introduction 

RADARSAT-2,  qui  devrait  etre  lance  vers  la  fin  de  2007,  transportera  un  radar  a  synthese 
d’ouverture  (RSO)  a  bande  C  qui  offrira  de  nombreux  modes  de  fonctionnement,  comme  celui  du 
RSO  polarimetrique  (PolSAR).  Un  systeme  PolSAR  emet  alternativement  des  impulsions 
electromagnetiques  a  polarisations  horizontale  (H)  et  verticale  (V),  puis  mesure  les  champs 
diffuses  a  polarisations  horizontale  (H)  et  verticale  (V).  II  existe  quatre  combinaisons  possibles  de 
champs  electriques  incidents  et  diffuses,  HH,  HV,  VH  et  VV.  RADARSAT-2  offrira  la 
quadrature  polarisation  (quad-pol),  la  polarisation  selective  (HH  et  HV  ou  VH  et  VV),  ainsi  que 
la  polarisation  simple  selective  (HH  ou  HV  ou  VH  ou  VV). 

Des  etudes  anterieures  ont  demontre  que  Tutilisation  des  donnees  de  PolSAR  peut  ameliorer 
considerablement  la  sensibilite  de  detection  des  navires  et  reduire  le  taux  de  fausses  alarmes, 
surtout  pour  les  petits  navires.  Les  donnees  a  quadruple  polarisation  procurent  la  plus  grande 
amelioration.  La  largeur  de  balayage  couverte  est  toutefois  reduite  en  raison  des  contraintes  de 
largeur  de  bande  du  systeme.  La  double  copolarisation  avec  phase,  qui  augmente  en  principe  la 
largeur  de  balayage  par  la  reduction  de  la  largeur  de  bande  requise,  permet  aussi  d’obtenir  un  bon 
rendement.  Le  rendement  se  caracterise  par  la  probability  de  detection  manquee  en  fonction  de  la 
probability  de  fausse  alarme. 

La  presente  etude  examine  la  possibility  d’appliquer  des  algorithmes  pour  ameliorer  la  detection 
et  la  classification  des  navires  par  Tutilisation  de  donnees  d’un  PolSAR  aeroporte  converties  en 
produits  de  polarimetrie  simules  de  RADARSAT-2.  L’algorithme  de  detection  des  navires  de 
PolSAR,  qui  a  ete  utilise  pour  cette  etude,  a  ete  mis  au  point  a  RDDC  Ottawa.  Trois  methodes  de 
classification  des  cibles  ont  ete  appliquees  :  les  methodes  de  Pauli,  de  Cameron  et  de 
caracterisation  de  diffusion  symetrique  (SSCM).  Le  logiciel  utilise  pour  les  methodes  de  Pauli  et 
de  Cameron  a  ete  mis  au  point  a  RDDC  Ottawa,  tandis  que  l’outil  SSCM  a  ete  incorpore  a  un 
poste  de  travail  PolSAR  developpe  au  Centre  canadien  de  teledetection  (CCT)  et  mis  a  la 
disposition  de  RDDC  Ottawa  par  l’Agence  spatiale  canadienne  (ASC). 

Les  images  de  RADARSAT-2  ont  ete  simulees  a  l’aide  du  simulateur  RADARSAT-2 
(RSAT2SIMU)  mis  au  point  par  MDA  et  ont  egalement  ete  mises  a  la  disposition  de  RDDC 
Ottawa  par  l’ASC.  Le  logiciel  RSAT2SIMU  simule  les  produits  RADARSAT-2  provenant  du 
RSO  polarimetrique  a  bande  C  CV- 580  d’Environnement  Canada  (EC),  dans  un  format  singulier 
complexe  (SLC),  en  elevant  correctement  le  seuil  de  bruit  et  en  degradant  la  resolution  spatiale. 

Le  simulateur  ne  peut  toutefois  pas  tenir  compte  de  la  geometric  particuliere  de  RADARSAT-2, 
c’est-a-dire  de  Taltitude,  de  Tangle  d’incidence  et  de  la  couverture  du  balayage.  La  qualite  des 
donnees  simulees  depend  de  la  qualite  des  donnees  de  RSO  CV- 580  fournies  a  T  entree. 
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Les  donnees  d’ entree  de  CV-580  utilisees  pour  ces  travaux  ont  ete  tirees  de  1’ experience  sur 
fintegration  des  capteurs  maritimes  (MARSIE).  L’ensemble  de  donnees  a  trait  a  plusieurs  navires 
connus,  notamment  le  Dominion  Victory  d’une  longueur  de  25  m. 

Portee 

L’etude  demontre  que  les  algorithmes  mis  au  point  pour  ameliorer  la  detection  et  la  classification 
des  navires  a  l’aide  des  donnees  de  PolSAR  aeroporte  s’appliqueront  aussi  aux  donnees 
polarimetriques  de  RADARSAT-2.  Les  resultats  mettent  en  evidence  1’ amelioration  significative 
du  rendement  de  detection  des  navires,  que  devrait  procurer  l’utilisation  de  donnees 
polarimetriques.  De  plus,  il  est  demontre  que  la  signature  polarimetrique  d’une  cible  peut 
ameliorer  le  rendement  de  classification  des  navires. 

Resultats 

Pour  la  detection  des  navires,  les  resultats  font  ressortir  l’avantage  des  systemes  polarimetriques. 
Par  exemple,  la  double  copolarisation  avec  phase  offre  un  meilleur  rendement  de  detection  que  la 
polarisation  HH  monocanal.  En  principe,  un  systeme  a  double  copolarisation  pourrait  produire 
une  plus  large  couverture  de  balayage  qu’un  systeme  polarimetrique. 

Pour  la  classification  des  cibles,  les  methodes  de  decomposition  des  cibles  de  Cameron  et  SSCM 
foumissent  de  1’ information  unique  sur  la  repartition  et  les  types  de  diffuseurs  qui  constituent  les 
cibles  d’interet,  offrant  ainsi  une  technique  de  classification  des  cibles.  Ces  methodes  de 
decomposition  s’averent  particulierement  efficaces  pour  les  gros  navires  bien  cibles.  La  methode 
de  decomposition  de  Pauli  offre  un  moyen  simple  mais  efficace  de  segmenter  la  cible  d’interet 
par  rapport  au  clutter  oceanique  environnant. 

Recherches  futures 

Les  resultats  demontrent  que  la  decomposition  des  cibles  et  l’analyse  des  signatures 
polarimetriques  peuvent  ameliorer  la  classification  des  cibles.  Pour  la  classification  des  cibles,  les 
aspects  suivants  sont  recommandes  en  vue  de  recherches  futures  : 

•  Relations  entre  les  proprietes  physiques  des  cibles  et  des  diffuseurs  elementaires,  selon  les 
methodes  de  decomposition  polarimetrique  de  Cameron  et  SSCM. 

•  Possibility  d’application  des  methodes  de  decomposition  polarimetrique  aux  algorithmes  de 
reconnaissance  automatique  des  cibles  (ATR)  par  l’utilisation  des  elements  de  diffusion  ou 
des  angles  estimes  d’ orientation  des  diffuseurs  comme  caracteristiques  des  cibles  dans  le 
logiciel  de  classification  des  cibles. 

•  Etablissement  d’algorithmes  fusionnant  toute  l’information  obtenue  des  methodes  d’analyse 
des  signatures  polarimetriques  pour  faciliter  la  classification  des  navires. 
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1.  Introduction 


Spacebome  C-band  polarimetric  synthetic  aperture  radar  capabilities  will  become  operationally 
available  in  the  near  future.  RADARS AT-2  is  expected  to  be  launched  in  2007  and  will  provide 
many  features,  including  quadrature  polarization  (quad-pol),  selective  polarization  (HH  and  HV 
or  VH  and  VV),  and  selective  single  polarization  (HH  or  HV  or  VH  or  VV)  [1],  where  the  H  and 
V  represent  horizontal  and  vertical  polarization  with  respect  to  the  antenna,  and  the  combinations 
represent  the  polarization  on  transmit  and  receive. 

In  order  to  support  the  development  of  PolSAR  ship  detection  and  classification  capabilities  for 
spaceborne  PolSAR  systems,  simulated  RADARS  AT-2  products  were  used  in  this  study  to 
investigate  the  feasibility  of  improved  ship  detection  and  classification  using  polarimetric  data. 
The  RADARS  AT-2  simulator  (RSAT2SIMU)  [2]  simulates  RADARS  AT-2  products  using  data 
from  the  Environment  Canada  (EC)  CV-580  C-band  PolSAR  in  single  look  complex  (SLC)  form 
by  increasing  the  noise  floor  and  degrading  the  spatial  resolution.  The  RSAT2SIMU  simulator 
was  developed  by  MDA  and  was  made  available  to  DRDC  Ottawa  by  CSA.  However,  the 
simulator  cannot  simulate  the  geometry  of  the  RADARS  AT-2,  such  as  the  altitude,  incidence 
angle  and  swath  coverage.  The  simulated  data  quality  depends  on  the  quality  of  the  input  CV-580 
data. 

In  this  report,  the  polarimetric  signatures  of  ships  from  simulated  Fine  Quad-Pol  and  Standard 
Quad-Pol  Modes  were  studied.  The  input  CV-580  data  was  obtained  during  the  MARSIE  trial 
[3]. 

1.1  Polarimetric  data  from  the  MARSIE  trial 

MARSIE  trial  is  a  “The  Technical  Cooperation  Program  (TTCP)”-lead  activity  that  is  exploring 
the  benefits  of  sensor  fusion  to  solve  the  target  detection  and  tracking  problem.  The  MARSIE 
trial  was  conducted  in  October  2005  off  the  East  Coast  of  Canada  and  brought  many  sensors  to 
bear  on  a  set  of  known  ship  targets  that  were  engaged  in  simulating  a  maritime  incursion  scenario 
(MIS). 

The  EC  CV-580  polarimetric  SAR  was  used  as  a  proxy  sensor  for  RADARS  AT-2  polarimetry,  as 
RADARS  AT-2  was  unavailable  during  the  trial  due  to  delays  in  its  launch.  MARSIE  provided  a 
unique  opportunity  to  collect  polarimetric  SAR  data  for  a  variety  of  known  ship  types  and  sizes 
carrying  out  realistic  manoeuvres. 

In  the  recent  report  on  the  processing  and  analysis  of  polarimetric  ship  signatures  from  MARSIE 
trial  for  project  Polar  Epsilon  [3],  the  results  clearly  demonstrated  the  significant  improvement 
that  may  be  realized  by  using  polarimetric  SAR  for  ship  detection  and  potentially  for  ship 
classification. 

This  report  provides  a  continuation  of  the  MARSIE  data  analysis,  focusing  mainly  on  simulated 
RADARS  AT-2  data  by  further  analyzing  the  ships  that  were  considered  in  [3]. 
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1.2  Methodology 


The  polarimetric  SAR  data  analysis  includes  target  detection  and  target  classification  assessment, 
similar  to  those  in  [3].  The  PolSAR  ship  detection  algorithm  used  in  this  study  was  developed  at 
DRDC  Ottawa  [4].  Three  polarimetric  decomposition  methods  were  used:  Pauli  [5],  Cameron  [6] 
and  SSCM  [7].  The  software  for  the  Pauli  and  the  Cameron  methods  was  developed  at  DRDC 
Ottawa,  while  the  SSCM  software  was  implemented  in  a  PolSAR  workstation  [9]  that  was 
developed  at  the  CCRS  and  was  made  available  to  DRDC  Ottawa  by  the  CSA. 

The  main  target  used  in  the  polarimetric  analysis  is  the  ship,  Dominion  Victory ,  which  was 
imaged  many  times  during  the  MARSIE  trial.  The  Dominion  Victory  has  a  length  of  25  m  which 
matches  the  minimum  detection  criteria  of  Project  Polar  Epsilon.  Two  other  ships  from  the 
MARSIE  trial,  Gulf  Service  and  HMCS  Toronto ,  were  also  studied.  The  lengths  of  these  ships  are 
42  m  and  134  m,  respectively. 

The  detection  performance  is  characterized  by  the  receiver  operating  characteristic  (ROC)  which 
shows  the  probability  of  missed  detection  as  a  function  of  the  probability  of  false  alarm.  The 
ROC  analysis  was  performed  for  various  radar  systems,  quadrature  polarization,  dual  co¬ 
polarization  (i.e.,  HH  and  VV)  with  both  amplitude  and  phase,  and  HH  single  polarization. 

The  classification  assessment  consists  of  applying  decomposition  methods  to  the  target  images  to 
determine  which  elemental  scatterers  from  each  method  are  present.  Histograms  of  the 
distributions  of  elemental  scatterers  were  studied  for  each  decomposition  method  and  each  target. 

For  the  Cameron  and  SSCM  methods,  the  results  were  compared  since  both  methods  determine 
the  maximum  number  of  symmetric  scatterers,  such  as  trihedrals  and  dihedrals. 

1 .3  Outline  of  this  document 

The  RADARS AT-2  capabilities  and  the  RADARS AT-2  simulator  are  described  in  Sections  2 
and  3,  respectively.  A  brief  overview  of  analysis  methods  is  presented  in  Section  4.  This  is 
followed  by  a  description  of  the  polarimetric  SAR  data  in  Section  5  and  the  data  analysis  results 
in  Section  6.  A  summary  of  the  results,  along  with  recommendations,  is  presented  in  Section  7. 
Many  details  have  been  relegated  to  Annexes  including  the  simulated  RADARS  AT-2  ship  images 
(Annex  A),  a  catalogue  of  polarimetric  decomposition  results  (Annexes  B  through  E),  and 
supporting  photographic  data  (Annex  F). 
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2.  RADARS AT-2  Capabilities 


The  RADARS  AT-2  modes  of  operation  and  swath  coverage  are  illustrated  in  Figure  1  [1]. 

The  Fine  Beam  Modes  are  intended  for  applications  that  require  higher  spatial  resolution  and  that 
can  be  successful  with  narrower  swath  coverage.  The  Fine  Beam  Modes  cover  an  incidence 
angle  range  from  30°  to  41°.  The  incidence  angle  is  the  angle  between  the  line  of  sight  from  the 
radar  and  the  local  vertical  direction. 

The  Standard  Beam  Modes  allow  imaging  over  a  wider  range  of  incidence  angles,  from  20°  to 
49°.  The  Sandard  Beam  Modes  offer  image  quality  characteristics  that  provide  a  compromise 
between  spatial  resolution,  radiometric  resolution,  and  swath  coverage. 

In  this  report,  we  focused  on  the  Fine  Quad-Pol  Mode  and  the  Standard  Quad-Pol  Mode.  The  key 
parameters  of  these  two  modes  are  summarized  in  Table  1. 


RADARSAT-2 


All  beam  modes 
available  in 
Left-looking  Mode 


Subsatellite  Track 

Beams 
(Low  Incidence) 


250 


Ultra- fine 
Beams 


Standard  Resolution 
Quad  Polarization 
(reduced  swath  width) 


Fine  Resolution 
Quad  Polarization 
(reduced  swath  width) 


Standard 

Beams 


Wide- 

Swath 

Beams 


Fine-Resolution 
Beams  (50  km  Swath) 


Multi-Look  Fine  Resolution 


Figure  1:  RADARSAT-2  operation  modes  and  coverage  (used  with  permission  [1]). 


Table  1:  Key  parameters  of  RADARSAT-2  polarimetric  image  modes. 


Beam  Mode 

Nominal 
Incidence 
Angle,  0inc  (°) 

Nominal 
Swath  Width 
(km) 

Nominal 

Swath 

Coverage 

(km) 

Approximate 
resolution 
GNDRg*  x  Az  (m) 

Quad-Pol 

(HH,  HV,  VH,  VV) 

Fine 

30-41 

25-50 

400  -  600 

11x9 

Standard 

20-41 

25-50 

250  -  600 

25  x28 

*  GNDRg:  Ground  range,  AZ:  Azimuth. 


DRDC  Ottawa  TM  2007-189 


3 


3.  RADARSAT-2  simulator 


The  RADARSAT-2  Simulator,  RSAT2SIMU,  is  intended  to  simulate  images  with  the  same 
resolution,  pixel  spacing  and  the  noise  equivalent  sigma  zero  (NESZ)  as  RADARSAT-2  products 
starting  from  EC  CV-5  80  PolSAR  data  [2], 

3.1  Simulation  process 

The  input  data  to  the  RSAT2SIMU  is  radiometrically  calibrated  EC  CV-580  SLC  data.  The 
simulation  of  the  RADARSAT-2  data  from  the  CV-580  data  involves  three  steps: 

•  Adding  zero-mean  complex  Gaussian  noise  to  increase  NESZ  to  RADARSAT-2  levels; 

•  Degrading  the  spatial  resolution  by  filtering  in  the  slant  range  and  azimuth  directions  to 
achieve  the  required  resolution  and  spectral  shapes,  while  preserving  polarimetric 
information;  and 

•  Re-sampling  by  interpolating  in  the  slant-range  and  azimuth  directions  to  achieve  the 
required  sample  spacing. 

The  key  system  parameters  of  the  CV-580  PolSAR  sensor  and  RADARSAT-2  are  compared  for 
single  look  data  in  Table  2.  The  range  sample  spacing  of  RADARSAT-2  is  calculated  using 
different  sampling  rates  for  various  modes,  and  in  azimuth  is  5  m,  which  is  about  10  times  larger 
than  that  of  the  CV-580.  The  range  resolution  is  in  slant  range  for  the  CV-580  and  in  ground 
range  for  RADARSAT-2.  The  slant  range  sample  spacing  listed  in  the  table  was  obtained  from 
the  output  header  files  of  the  simulated  RADARSAT-2  images.  Each  polarimetric  channel  was 
processed  separately  using  the  same  parameters.  All  four  channels  (HH,  HV,  VH,  and  VV)  were 
processed  for  the  polarimetric  modes. 

Table  2:  Key  system  parameters  for  the  CV-580  and  RADARSAT-2  SARs. 


CV-580 

RADARSAT-2 

Quad  Fine 

Quad  Standard 

Frequency  (GHz) 

5.3 

5.405 

Altitude  (km) 

5  to  7 

798 

Azimuth  resolution  (m) 

~1 

9 

28 

Range  resolution  (m) 

Slant-range 

6 

Ground-range 

11 

25 

Azimuth  sample  spacing  (m) 

0.4 

5.1 

Slant  range  sample  spacing  (m) 

4 

4.7* 

11.9* 

NESZ  (dB) 

~  -40 

-31 

-28 

Platform  speed  (m/s) 

~140 

7460 

*  According  to  the  output  header  file  of  the  simulated  image. 


RSAT2SIMU  provides  both  a  RADARSAT-2  product  format  and  an  intermediate  product  in  CV- 
580  data  format.  In  this  report,  the  simulated  RADARSAT-2  images  in  the  CV-580  data  format 
were  used. 
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As  an  example,  an  input  CV- 580  image  and  the  simulated  RADARSAT-2  Fine  Quad-Pol  and 
Standard  Quad-Pol  Mode  images  in  CV-580  format  are  illustrated  in  Figure  2  and  Figure  3, 
repectively.  This  image  includes  an  urban  environment,  residential  areas,  buildings,  parks, 
bridges,  bodies  of  water,  etc.  Figure  3a)  and  3b)  are  intended  to  be  illustrative  of  the  simulated 
RADARSAT-2  Fine  Quad-Pol  and  Standard  Quad-Pol  images  from  the  CV-580  data. 


>  Azimuth 


▼ 

Range  Figure  2:  CV-580  data  input  to  RSAT2SIMU. 


Figure  3:  Simulated  RADARSAT-2  images;  a)  Fine  Quad-Pol  Mode:  b)  Standard  Quad-Pol  Mode. 


3.2  Differences  with  respect  to  RADARSAT-2  products 

There  are  a  number  of  differences  between  the  simulated  products  and  actual  RADARSAT-2 
products  relevant  to  this  study: 

•  The  polarimetric  SAR  data  quality  has  significant  effects  on  the  polarimetric  signatures  of 
the  targets  and  the  data  quality  depends  largely  on  the  various  calibrations,  such  as 
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radiometric,  phase  and  geometric  calibration  [17].  The  simulated  data  quality  is  determined 
by  the  CV-580  calibration  system  and  not  the  RADARSAT-2  calibration  system; 

•  The  NESZ  is  assumed  constant  over  both  the  input  and  output  images.  In  practice,  NESZ 
varies  over  the  swath,  being  lower  at  the  centre  of  the  swath  and  higher  at  the  swath  edges. 
This  variation  is  caused  by  the  antenna  pattern  of  the  sensor,  which  varies  from  beam  to 
beam;  and 

•  The  simulated  products  do  not  represent  a  particular  RADARSAT-2  beam.  The  CV-580 
covers  a  wide  range  of  incidence  angles,  far  exceeding  the  incidence  angle  range  of  the 
RADARSAT-2  beams.  CV-580  coverage  is  limited  compared  to  that  of  RADARSAT-2. 
The  size  of  the  simulated  region  does  not  represent  the  size  of  a  RADARSAT-2  product. 
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4.  Overview  of  analysis  methods 


The  polarimetric  SAR  data  analysis  includes  ship  detection  and  target  classification  assessment. 
For  ship  detection,  we  studied  the  detection  performance  of  various  PolSAR  systems  and 
estimated  the  receiver  operating  characteristic  (ROC).  The  performance  of  three  PolSAR  system 
configurations  was  compared:  polarimetric,  dual  co-polarization  (i.e.,  HH  and  VV)  with  relative 
phase,  and  single  polarization  HH. 

For  target  classification,  three  polarimetric  decomposition  methods  were  considered:  Pauli, 
Cameron,  and  SSCM  ([5],  [6],  [7]). 

4.1  Receiver  operating  characteristic 

The  ROC  is  plotted  as  the  probability  of  missed  detection  {PMd)  versus  probability  of  false  alarm 
(Pfa)-  Measured  data  processed  by  COASP  [10]  were  used  to  calculate  these  probabilities  using  a 
well-known  statistics-based  methodology  based  upon  likelihood  ratio  tests  with  the  Neyman- 
Pearson  criterion  ([4],[1 1]). 

A  polarimetric  SAR  system  provides  observation  of  the  scattering  matrix  X(i,  j )  for  each  pixel 
(i,  j )  in  an  image.  The  matrix  components  SHH,  SHV,  SVH ,  and  Sw  are  complex  valued  elements, 
representing  the  scattering  magnitude  and  phase  for  the  possible  combinations  of  horizontal  (H) 
or  vertical  (V)  transmit  and  receive  polarizations.  The  components  of X(i,  j )  can  be  written  as  a 
vector: 

X(i,j)  =  [SHH(i'>j)  SHV(iJ)  SVH  (/,  j)  Syy(i,j)]T  (1) 


where  the  superscript  T  is  the  transpose  operator. 

Ship  detection  is  a  binary  decision  problem.  The  fundamental  algorithms  of  PolSAR  ship 
detection  are  based  on  statistical  decision  theory  that  may  be  applied  directly  to  the  components 
of  the  scattering  matrix  to  obtain  a  decision  variable.  A  likelihood  ratio  test  with  the  Neyman- 
Pearson  criterion  is  used  to  define  a  point  detection  criterion.  The  Gaussian  distributions  for 
scattering  matrix  components  are  assumed  to  derive  an  approximate  decision  variable,  while  the 
measured  data  may  be  used  to  calculate  the  detection  variables. 


Following  [4],  the  decision  variable  is  approximately  given  by 


Xi/(C;1)X  = 


>  TJ 


[<ri 


for  a  ship 


for  ocean 


(2) 


where  C0  is  the  covariance  matrix  of  the  ocean. 

For  a  polarimetric  system,  X  is  given  by  Equation  (1);  for  a  dual  co-polarization  system,  such  as 
HH/VV  with  amplitude  and  phase,  X  is  given  by 
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(3) 


X(z,y)  -[ShH&J)  Syy(i9j)]T 


For  a  single  polarisation  system,  the  decision  variable  is  given  by  the  amplitude  as 


XI 


\ 

r 

1 

_  < 

2 

) 

>ri 

<?l 


for  a  ship 


for  ocean 


(4) 


where  <j0  is  the  mean  value  of  the  ocean  intensity.  From  these  algorithms,  we  can  see  that  a 
polarimetric  system  contains  amplitude  and  phase  information  from  four  channels.  As  such,  it  is 
expected  to  provide  the  best  detection  performance  compared  to  systems  with  fewer  channels.  A 
dual  co-polarization  system  with  amplitude  and  phase  contains  information  from  two  channels,  so 
it  should  provide  better  detection  performance  than  a  single  channel  system  which  contains  only 
the  amplitude  information. 

The  detection  performance  is  determined  by  two  factors:  the  probability  of  false  alarm  (PFA)  and 
the  probability  of  missed  detection  (Pmd)-  The  ROC  plots  the  probability  of  missed  detection 
versus  the  probability  of  false  alarm.  Measured  data  is  used  to  calculate  these  probabilities  as  the 
detection  threshold  rj  is  varied. 

The  algorithms  for  calculation  of  the  ROC  can  be  applied  to  any  designed  probability  of  false 
alarm,  such  as  PFA  =  10'9,  provided  there  are  enough  ocean  samples  available.  The  algorithms 
are  independent  of  the  environmental  conditions;  they  can  be  applied  to  any  sea  state  condition. 

The  accuracy  of  the  estimated  values  of  PMd  and  PFA  may  be  determined  from  the  number  of 
missed  detection  and  false  alarm  events.  For  a  PFA  =  10"9,  to  obtain  an  error  of  about  30%  in  the 
estimate  of  PFA,  roughly  10 10  independent  ocean  samples  at  the  same  geometry  and  environmental 
conditions  are  required;  to  obtain  an  error  of  about  10%,  roughly  1011  independent  ocean  samples 
are  required.  Since  these  numbers  are  rather  large,  PFA  =  10"4  is  used  to  estimate  the  relative 
performance,  which  for  106  ocean  samples  gives  an  error  of  10%. 

For  target  detection,  several  steps  are  required.  In  this  study,  the  detection  is  performed  on  each 
pixel  in  the  image  and  PFA  =  10"4  is  set  for  considerations  of: 

•  Obtaining  lower  PMD\  and 

•  Constraining  the  estimation  accuracy  due  to  the  limited  number  of  available  samples. 

There  is  a  trade  off  between  PFA  and  PMd •  In  general,  the  lower  PFA  is  set  for  detection,  the  higher 
Pmd  will  be.  To  obtain  a  lower  PFA  with  a  lower  PMd ,  a  combination  of  pixel  detection  with  other 
algorithms  must  be  employed.  For  example,  following  pixel-based  detection,  other  algorithms 
such  as  target  clustering,  sub-aperture  analysis,  or  polarimetric  signature  analysis  may  be  applied 
to  achieve  a  lower  PFA  with  a  low  PMd -  For  clustering,  the  detected  pixels  must  be  grouped  and  a 
decision  made  as  to  whether  or  not  there  are  sufficient  pixels  in  the  group  to  represent  a  target. 
Such  a  procedure  will  improve  PFA  and  PMd ,  but  at  the  cost  of  increasing  the  minimum  size  of 
target  that  can  be  detected.  Phase  interferogram  analysis  between  the  HV  and  VH  channels  in 
polarimetric  systems  can  aid  in  distinguishing  the  targets  from  the  ambiguities  [12].  This  will 
also  reduce  the  PFA. 
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4.2 


Polarimetric  target  decomposition 


4.2.1  Pauli 


A  linear  polarization  scattering  matrix  S  can  be  expressed  in  the  Pauli  basis  [5]  as: 


V2 


1  0 
0  1 


s 


1  0 
0  -1 


0  1 
1  0 


V2 


0  -i 
i  0 


(5) 


This  representation  permits  the  extraction  of  physical  information  from  the  2x2  coherent 
scattering  matrix.  Under  the  reciprocity  assumption  (i.e.,  SHv  =  Svh),  the  Pauli  scattering  vector  is 
given  by  [5]: 


vv 


$HH ~  Syv 


Shv+Svh] 


(6) 


In  Pauli  decomposition  analysis,  each  pixel  is  assigned  to  one  of  three  classes,  depending  on 
which  of  the  components  is  the  largest  in  amplitude.  For  class  1,  (SHH  +  Svv )  is  the  largest, 
representing  single  or  odd  bounce  scattering;  for  class  2,  (.S'////  Sn  )  is  the  largest  representing 
double  or  even  bounce  scattering;  and  for  class  3,  ( SHv  +  Svh)  is  the  largest  representing  volume 
scattering. 

The  physical  interpretation  of  three  basic  scattering  mechanisms  is  illustrated  in  Figure  4. 


SINGLE  BOUNCE  DOUBLE  BOUNCE 

SCATTERING  SCATTERING 


VOLUME 

V.  SCATTERING 


HV  +  VH 


Figure  4:  Physical  interpretation  of  three  basic  scattering  mechanisms. 


4.2.2  Cameron 

In  the  Cameron  method  [6],  the  measured  scattering  matrix  is  decomposed  based  on  two  basic 
properties  of  radar  returns:  reciprocity  and  symmetry.  An  arbitrary  scattering  matrix  could  be 
coherently  decomposed  into  a  nonreciprocal  component,  a  maximum  symmetric  component,  and 
a  minimum  symmetric  component.  In  this  report,  only  the  maximum  symmetric  scattering 
components  are  studied.  These  components  are  trihedral,  diplane  (i.e.,  dihedral),  dipole,  cylinder, 
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narrow  diplane,  and  quarter-wave  devices.  The  symmetric  scatterers  are  represented  on  a  Unit 
disk  as  shown  in  Figure  5. 


Figure  5:  Unit  disk  representation  of  symmetric  scatterers. 


A  symmetric  scatterer  can  be  represented  in  terms  of  its  normalized  diagonal  scattering  matrix 

A(z): 


A(z)  = 


rn 

0 

0 


VZ7 


Izl  <  1 


(7) 


where  the  parameter  z  determines  the  scatterer  type  associated  with  a  particular  scattering  matrix, 
which  can  be  represented  by  a  point  in  the  unit  disc  of  the  complex  plane.  The  six  elemental 
scatterers  are  expressed  in  terms  of  A(z)  as  follows: 


trihedral 

diplane/ 

dihedral 

dipole 

cylinder 

narrow 

diplane 

quarter 

wave 

A(l) 

A(-l) 

A(0) 

A(‘/2) 

M-K) 

A(±i) 

Each  pixel  is  then  assigned  to  one  of  six  classes,  depending  on  which  of  the  components  is  the 
shortest  distance  d  away,  according  to  the  symmetric  scatter  distance  measure: 


d(zi,Z2)  =  cos 


maxQl  +  z^  \,\zl+zl\) 

I  Z1  I2  |  z2  I2 


0<d  <  — 

K 


(8) 
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where  z\  and  z2  are  the  z  values  for  the  two  scatterers  under  consideration.  In  general,  trihedral 
represents  odd  bounce  scattering  and  diplane  represents  even  bounce  scattering,  as  in  the  Pauli 
decomposition.  Also,  cylinder  represents  curved  metallic  plates  as  well  as  large  diameter  or  pipe¬ 
like  objects. 

4.2.3  Symmetric  Scattering  Characterization  Method 

The  SSCM  is  another  approach  to  exploit  the  information  from  the  largest  target  symmetric 
scattering  components  [7].  These  components  are  the  same  as  those  defined  in  the  Cameron 
method. 

The  symmetric  scattering  vectors  in  SSCM  are  represented  in  terms  of  the  latitude  2\|/  and 
longitude  2%  on  a  target  Poincare  sphere  surface,  as  shown  in  Figure  6,  instead  of  on  the  Unit 
Disk  considered  in  the  Cameron  method.  The  coordinates  of  the  six  largest  target  symmetric 
elemental  scatterers  on  the  Poincare  sphere  are  summarized  in  Table  3. 


1/4  Wave 


Figure  6:  Poincare  sphere  representation  of  symmetric  scatterers  (use 

with  permission  [7]). 


Table  3.  Elemental  scatterer  coordinates  on  the  Poincare  sphere. 


scattering 

trihedral 

dihedral 

dipole 

cylinder 

narrow 

diplane 

quarter 

wave 

(VoZc) 

(0°,  0°) 

(90°,  0°) 

(45°,  0°) 

(18.43°,  0°) 

(81.88°,  0°) 

(0°,  45°) 

Based  on  the  symmetric  scatterer  centre  on  the  Poincare  sphere,  six  symmetric  scatterers  can  also 
be  generated  by  using  the  thresholds  shown  in  Table  4. 
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For  example,  a  symmetric  scatterer  is  classified  as  a  trihedral  scatterer  if  its  Poincare  sphere 
longitude  angle  (//.  is  between  0°  and  9.22°  and  its  Poincare  sphere  latitude  angle  %c  *s  between 
-12.5°  and  +12.5°. 


Table  4:  threshold  to  generate  each  symmetric  scatterer. 


symmetric  scatterer 

Vci°) 

XcO 

trihedral 

0  to  9.22 

and 

-12.5  to  12.5 

cylinder 

9.22  to  31.72 

dipole 

31.72  to  63.44 

narrow  diplane 

63.44  to  85.94 

dihedral 

85.94  to  90 

quarter  wave 

Oto  90 

and 

<-12.5  or  >  12.5 
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5  Polarimetric  SAR  data 


5.1  Data  description 

The  PolSAR  data  used  in  this  report  are  listed  in  Table  5.  They  were  obtained  during  the 
MARSIE  [3]  and  the  Urban  PolInSAR  [14]  trials  using  the  EC  CV-580  PolSAR  sensor. 

Table  5:  CV-580  data  used  in  the  RADARSAT-2  simulations. 


Flight  line/pass,  date 

Sensor  look  direction 

Trial  name 

llp8,  04-Sep-02 

left 

Urban  PolInSAR 

141pl,  20-0ct-05 

Right 

MARSIE 

141p2, 20-0ct-05 

Left 

MARSIE 

141p3, 20-0ct-05 

Right 

MARSIE 

142p4,  20-0ct-05 

left 

MARSIE 

122p2,  17-Oct-05 

left 

MARSIE 

The  Urban  PolInSAR  trial  data  were  processed  using  the  Configurable  Airborne  SAR  Processor 
(COASP)  since  the  targets  of  interest  in  this  image  are  static  targets.  However,  the  individual 
ship  images  from  the  MARSIE  trial  were  reprocessed  using  the  Chip-based  Adaptive  SAR 
Processor  (CHASP)  to  improve  the  ship  focus  by  compensating  for  the  ship  motion.  Both 
processors,  CO  ASP  and  CHASP,  were  developed  at  DRDC  Ottawa  [10]. 

The  Urban  PolInSAR  data  were  used  to  evaluate  the  quality  of  the  simulated  data,  while  the  main 
targets  used  in  the  polarimetric  analysis  are  vessels  which  were  obtained  during  the  MARSIE 
trial,  as  listed  in  Table  6.  The  length  of  the  Dominion  Victory  matches  the  minimum  detection 
criteria  of  project  Polar  Epsilon;  it  was  imaged  several  times  from  various  incidence  and  aspect 
angles.  The  aspect  angle  is  the  angle  from  the  orientation  of  the  ship  to  the  Azimuth  direction 
(clockwise). 

Table  6:  Ships  used  in  the  analysis 


Ship 

Length  (m) 

Type 

Dominion  Victory 

25 

Multi-purpose  support 

Gulf  Service 

42 

Tug 

HMCS  Toronto 

134 

Halifax  Class  Patrol  Frigate 

5.2  Data  quality  evaluation 

In  order  to  evaluate  the  quality  of  simulated  RADARSAT-2  data,  a  known  corner  reflector  (CR), 
referred  to  as  S5  and  deployed  at  the  Cannaught  Range  calibration  site  in  the  Ottawa  area,  was 
studied.  The  HH  images  of  the  calibration  site  from  the  CV-580  and  simulated  RADARSAT-2 
Fine  Quad-Pol  and  Standard  Quad-Pol  Modes  are  shown  in  Figure  7.  The  CV-580  image  was 
obtained  from  flight  line  1  pass  8  (llp8)  during  the  Urban  PolInSAR  trial.  Comer  reflector  S3 
has  a  70.5  cm  back  spine  length. 
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Since  the  HH  channel  image  from  RADARSAT-2  should  be  similar  to  a  RADARSAT-1  image,  a 
RADARSAT-1  image  obtained  from  7  May  2004  is  shown  in  Figure  8  as  a  comparison.  In 
addition,  a  corner  reflector  having  a  back  spine  length  of  75  cm,  located  at  N45. 37368  and 
W75. 92441  in  the  image,  was  analyzed;  the  results  were  compared  to  those  from  the  simulated 
data  since  these  CRs  have  similar  size.  This  image  is  a  Fine  2  Near  beam  mode  with  a  nominal 
resolution  of  9  m  in  ground  range  and  9  m  in  azimuth,  acquired  in  the  ascending  pass  direction 

[15]. 
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Figure  7:  Corner  reflectors  and  calibration  site  in  HH  images  (llp8,  24  September  2002);  a) 
CV-580;  b)  simulated  RADARSAT-2  Fine  Quad-Pol  Mode;  c)  simulated  RADARSAT-2 

Standard  Quad-Pol  Mode. 
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Figure  8:  Calibration  site  and  75  cm  corner  reflector  in  a  RADARSAT-1  image  (7  May  2004). 


The  analysis  results,  such  as  peak  value,  mean  clutter  value,  peak-to-clutter  ratio  (PCR),  and 
resolutions  derived  from  the  -3dB  “width”  response  of  the  CR,  for  the  various  modes  are 
summarized  in  Table  7.  The  target  peak  to  clutter  ratio  (PCR)  of  Standard  Quad-Pol  Mode 
simulated  RADARS AT-2  HH  data  is  almost  the  same  as  that  of  actual  RADARSAT-1  data,  but 
the  simulated  Fine  Quad-Pol  Mode  PCR  is  much  higher. 

The  simulated  RADARS  AT-2  product  co-pol  and  cross-pol  channels  are  also  compared.  The  HH 
and  VV  corner  reflector  measurements  results  are  similar  while  HV  and  VH  responses  are  much 
lower. 

The  resolutions  of  the  simulated  RADARS  AT-2  images  in  the  azimuth  direction  are  comparable 
to  the  actual  resolution  (see  Table  1).  It  is  worth  noting  that  the  azimuth  resolution  of  28  m  (see 
Table  1)  in  the  actual  RADARSAT-2  Standard  Quad-Pol  Mode  is  for  4-looks.  However,  the 
azimuth  resolution  of  the  simulated  products  is  for  one-look  (see  Table  7).  The  range  resolution 
of  the  simulated  images  can  not  be  compared  to  the  actual  products  since  the  range  represents  the 
slant  range  in  the  simulated  data,  but  represents  the  ground  range  in  the  actual  data. 
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Table  7:  Statistical  parameters  of  a  corner  reflector  and  the  image  resolution. 


S5 

Peak 

(dB) 

Mean  Clutter 
(dB) 

PCR 

(dB) 

Resolution  (m) 

CV-580 

HH 

16.40 

-17.99 

34.39 

6x1,  Slant  Range  x  AZ* 

Simulated 

RADARSAT-2 

FQ* 

HH 

6.54 

-18.48 

25.02 

5.88  x  8.92,  Slant  Range  x  AZ 

StdQ* 

HH 

3.23 

-17.57 

20.80 

14.63  x  9.59,  Slant  Range  x  AZ 

RADARSAT-1 

HH 

8.38 

-11.90 

20.28 

9X9,  Ground  Range  x  AZ 

Simulated 

RADARSAT-2 

FQ 

VV 

5.36 

-18.59 

23.95 

StdQ 

vv 

2.01 

-18.54 

20.56 

Simulated 

RADARSAT-2 

FQ 

HV 

-15.11 

-23.99 

8.88 

StdQ 

HV 

1.55 

-25.47 

27.01 

Simulated 

RADARSAT-2 

FQ 

VH 

-13.83 

-24.59 

10.76 

StdQ 

VH 

1.07 

-25.52 

26.59 

*  FQ:  Fine  Quad-Pol  Mode,  StdQ:  Standard  Quad-Pol  Mode,  AZ:  Azimuth. 


Furthermore,  an  urban  image,  consisting  of  buildings,  streets,  bridges,  parks,  etc.,  was  also 
studied.  The  images  from  the  HH  channel  of  simulated  RADARS AT-2  Fine  Quad-Pol  and 
Standard  Quad-Pol  Modes  and  RADARSAT-1  are  shown  in  Figure  9  and  Figure  10,  respectively. 
In  order  to  compare  these  images,  the  simulated  RADARS  AT-2  images  are  flipped  vertically  due 
to  the  different  acquisition  geometry. 

The  urban  structures  observed  in  the  simulated  RADARS  AT-2  images  are  similar  to  those  in  the 
RADARSAT-1  image.  However,  the  simulated  RADARSAT-2  images  will  differ  from  actual 
images,  as  discussed  in  Section  3.2. 

From  these  analysis  results,  the  simulated  RADARSAT-2  images  are  considered  to  be  sufficient 
for  the  evaluation  of  the  application  of  polarimetric  SAR  analysis  algorithms  developed  using 
airborne  PolSAR  data  to  spaceborne  systems,  specifically  for  ship  detection  and  target 
classification  methods. 
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%K 


b) 

Figure  9:  Simulated  RADARSAT-2  images  of  urban  environment  ( llp8 ,  24  September  2002); 

a)  Fine  Quad-Pol  Mode;  b)  Standard  Quad-Pol  Mode. 
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Figure  10:  RADARSAT-1  image  of  urban  environment  (7  May  2004). 
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6  Analysis  results 


The  polarimetric  SAR  data  analysis  includes  target  detection  and  target  classification  assessment. 
The  target  detection  performance  is  estimated  using  the  ROC  curve  and  the  classification  is 
carried  out  using  three  polarimetric  decomposition  methods  and  the  target  signature  analysis. 

For  ROC  analysis,  the  standard  PolSAR  processed  data  were  used.  However,  the  most  standard 
PolSAR  processed  ships  are  smeared  due  to  ship  motion.  Therefore,  for  target  classification 
analysis,  the  individual  ships  were  reprocessed  to  improve  the  focus  before  RADARS AT-2 
simulation. 

Examples  of  the  analysis  results  of  the  Dominion  Victory  from  141pl  are  presented  in  this 
Section;  more  results  can  be  found  in  Annexes  A  through  E. 

6.1  Simulated  vessel  images 

An  example  of  a  CV-580  CHASP-processed  image  and  the  simulated  RADARS  AT-2  Fine  Quad- 
Pol  Mode  and  Standard  Quad-Pol  Mode  images  are  shown  in  Figure  11.  The  ship  is  visible  in 
both  simulated  RADARSAT-2  images.  However,  in  141p2,  the  ship  doesn’t  appear  in  the 
simulated  images  (See  Annex  A). 


v  100  200  300  400  500  600  700  800  900  1000  1100 

a) 
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Figure  11:  Polarimetric  image  of  Dominion  Victory  (I41pl,  20  October  2005);  a)  CV-580;  b) 
simulated  RADARSAT-2  Fine  Quad-Pol  Mode;  c)  simulated  RADARSAT-2  Standard 

Quad-Pol  Mode. 


6.2  Receiver  Operating  Characteristic 

The  detection  performance  was  estimated  for  three  systems,  polarimetric,  dual  co-polarization 
with  amplitude  and  phase,  and  single  channel  HH.  The  results  are  shown  in  Figure  12. 

The  advantage  of  a  quad  polarimetric  system  is  apparent  (i.e.,  decreased  probability  of  missed 
detection,  decreased  probability  of  false  alarm),  which  uses  the  available  amplitude  and  phase 
information  (triangles),  in  both  Find  Quad-Pol  and  Standard  Quad-Pol  Modes.  The  dual  co¬ 
polarization  system  with  phase  (circles)  provides  better  detection  performance  than  a  single 
channel  HH  polarization  system  (asterisks).  In  principle,  a  dual  polarization  system  should 
provide  wider  swath  coverage  than  a  quad  polarimetric  system.  This  performance  improvement 
as  polarimetric  channels  are  added  has  been  demonstrated  previously  (see  [3],  [4]  and  Annex  B). 

The  variability  in  ship  detection  of  various  systems  is  also  show  in  Figure  12a.  For  this  example, 
the  Fine  Quad-Pol  Mode  had  an  ocean  sample  composed  of  630,000  pixels  and  a  ship  sample 
composed  of  160  pixels,  while  the  Standard  Quad-Pol  Mode  had  an  ocean  sample  of  227,000 
pixels  and  ship  sample  composed  of  96  samples.  The  available  number  of  ocean  samples 
depends  on  the  maximum  data  size  that  can  be  handled  by  the  simulator,  while  the  maximum 
number  of  ship  samples  depends  on  the  size  of  the  ship  and  the  sample  spacing  of  the 
RADARSAT-2  data. 

By  arbitrarily  selecting  a  probability  of  false  alarm  of  PFA  =  10“4,  the  relative  improvement  in  ship 
detectability  may  be  quantified  by  comparing  the  probability  of  missed  detection  PMD  across  the 
cases  studied.  The  relative  improvement  in  the  detection  performance  of  all  three  images  from  20 
October  2005  is  summarized  in  Table  8  by  taking  the  ratio  of  the  observed  PMd  to  that  of 
polarimetry  PQmd •  In  Table  8,  entries  where  no  number  is  given  indicates  no  detection.  From 
Table  8,  we  can  see  that  in  some  cases,  Dominion  Victory  could  not  be  detected  using  a  single 
channel  system.  Furthermore,  it  could  not  always  be  detected  even  when  using  a  dual  co-pol 
system.  The  probability  of  missed  detection  is  often  high  for  single  channel  and  dual-polarized 
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systems.  The  results  demonstrate  that  the  polarimetric  system  is  always  the  best.  The  dual  co¬ 
polarized  HH-VV  system  provides  better  detection  than  a  single  channel  HH  polarization  system. 


io'6  io'5  io'4  io'3  io'2  10'1  10° 

False  alarm  probability 


358  DV,  R-2  Simulations  Standard 


10'6  10'5  10'4  10'3  10'2  10'1  10° 
False  alarm  probability 


358  DV,  R-2  Simulations  Fine  Beam 


a  b 

Figure  12:  Detection  performance  for  Dominion  Victory  ( 141pl ,  20  October  2005).  a)  Fine 

Quad-Pol  Mode;  b)  Standard  Quad-Pol  Mode. 


Table  8:  Values  of  PMd  for  PFa  =  1  O'4  for  Dominion  Victory. 


Dominion  Victory  (L=25m) 

line/pass 

Find  Quad-Pol 

Standard  Quad-Pol 

PQmd 

Pmd/PQmd 

PQmd 

Pmd/PQmd 

14  lpl 
20-0ct-05 

Quad 

Pmd 

0.15 

1.00 

0.34 

1.00 

HH-VV 

Pmd 

0.22 

1.47 

0.45 

1.32 

HH 

Pmd 

0.28 

1.87 

0.49 

1.44 

14  lp2 
20-0ct-05 

Quad 

Pmd 

0.23 

1.00 

0.60 

1.00 

HH-VV 

Pmd 

- 

- 

0.70 

1.17 

HH 

Pmd 

- 

- 

- 

- 

142p4 

20-0ct-05 

Quad 

Pmd 

0.38 

1.00 

0.56 

1.00 

HH-VV 

Pmd 

0.49 

1.29 

- 

- 

HH 

Pmd 

0.70 

1.84 

- 

- 

6.3  Target  decomposition 

The  three  target  decomposition  methods,  Pauli,  Cameron  and  SSCM,  were  applied  to  the 
simulated  RADARSAT-2  data.  For  the  Pauli  and  Cameron  decomposition  analysis,  a  threshold 
for  each  ship  was  selected  based  on  the  estimated  mean  ocean  clutter  value  (a°)  to  reduce  the 
clutter  in  the  decomposition  results.  The  thresholds  are  listed  in  Table  9. 
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Table  9:  Decomposition  analysis  thresholds. 


Line/Pass 

Ship  Name 

Reference  Channel 

Threshold  (dB) 

Fine  Quad-Pol  and 
Standard  Quad-Pol 

17  Oct.  2005 

L22p2  (a342) 

Gulf  Service 

HH 

-15 

122p2  (a342) 

HMCS  Toronto 

HH 

-20 

20  Oct.  2005 

141pl  (a358) 

Dominion  Victory 

HH 

-15 

141p3  (a360) 

Dominion  Victory 

HH 

-12 

142p4  (a361) 

Dominion  Victory 

HH 

-14 

6.3.1  Pauli 


The  Pauli  decomposition  results  are  shown  in  Figure  13,  in  which  each  pixel  is  assigned  to  one  of 
three  classes,  depending  on  which  of  the  components  has  the  largest  amplitude.  For  Class  1,  (SHh 
+  *Syv)  is  the  largest,  for  class  2,  (SHh  -  *SW)  is  the  largest,  and  for  class  3,  (Shv  +  SVh)  is  the 
largest.  Class  1  represents  odd  bounce  (blue),  Class  2  represents  even  bounce  (red),  and  the  class 
3  represents  volume  scattering  (green).  As  might  be  expected,  the  ship  image  includes  mainly 
double  bounce  scattering  and  single  bounce  scattering  due  to  ship  structure.  Volume  scattering 
also  appears  in  the  ship  image  in  some  cases.  This  may  be  due  to  multiple  bounces  from  the  ship 
structure  or  may  be  caused  by  uncompensated  ship  motion  during  image  acquisition  and 
processing. 

It  is  noted  that  there  are  more  single  bounce  scattering  components  in  the  simulated  images  than 
in  the  original  CV-580  images  ([3]). 


a)  b) 

Figure  13:  Pauli  decomposition  images  of  Dominion  Victory  (I41pl  20  October  2005).  a)  Fine 

Quad-Pol  Mode;  b)  Standard  Quad-Pol  Mode. 
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6.3.2  Cameron 


The  Cameron  decomposition  results  are  shown  in  Figure  14.  The  ship  image  has  been 
decomposed  into  six  elemental  scatterers.  In  this  image,  the  most  dominant  component  is  the 
cylinder-like  scatter,  followed  by  the  dipole  and  the  quarter-wave  device.  The  distribution  of 
each  component  is  shown  in  the  histogram.  The  types  of  elemental  scatterers  are  indicated  as 
follows:  1:  none  classified,  2:  trihedral,  3:  diplane  (i.e.  dihedral),  4:  dipole,  5:  cylinder,  6:  narrow 
diplane  and  7:  quarter- wave  device. 

For  the  cases  studied,  the  dominant  component  varies.  However,  the  distribution  of  dominant 
components  from  the  Standard  Quad-Pol  Mode  is  similar  to  those  from  the  Fine  Quad-Pol  Mode. 
The  fractions  of  each  elemental  scatterer  from  the  Cameron  decomposition  are  summarized  in 
Table  10  and  Table  11.  This  report  has  focused  on  the  distributions  of  numbers  of  elemental 
scatterers  in  the  target  displayed  as  a  histogram.  The  distribution  of  scatterers  is  potentially 
useful  for  ship  classification.  But  it  is  recommended  to  combine  this  information  with  other 
target  analysis  methods.  It  is  judged  that  these  distributions  are  not  always  different  enough  to 
permit  clear  classifications.  Therefore,  the  use  of  the  spatial  distributions  of  scatterers  should  also 
be  investigated. 


358  Dominion  Ship 


358  Dominion  Ship 


11/4 
ndpln 
cyl 
dpi 
dpln 
-tri 

-  none 


358  Dominion  Ship  l.none,  2.tri,  3. dpln,  4.dpl,  5.cyl,  B.ndpln,  7.1/4 
0.35 
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0.2 
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1  2  3  4  5  6  7 


Figure  14:  Cameron  decomposition  image  (left)  and  histogram  (right)  of  Dominion  Victory 
(I41pl  20  October  2005).  Top:  Fine  Quad-Pol ,  Bottom:  Standard  Quad-Pol. 

In  the  histograms  (right):  the  scatterer  type  is  along  the  x-axis  and  the  percentage  of  the 

distribution  is  along  the  y-axis. 
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Table  10:  Summary  of  elemental  scatter  er  distributions  for  Cameron  decomposition  for 

simulated  Fine  Quad-Pol  images. 


Ship  name 

trihedral 

diplane 

dipole 

cylinder 

narrow 

diplane 

quarter 

wave 

17  Oct.  2005, 122p2  (a342)  1 

Gulf  Service 

0.1811 

0.0157 

0.0728 

0.2953 

0.0296 

0.4055 

HMCS  Toronto 

0.0386 

0.0289 

0.3269 

0.2251 

0.1715 

0.2090 

20  Oct.  2005,  Dominion  Victory 

141pl  (a358) 

0.1114 

0.0107 

0.2672 

0.3023 

0.0718 

0.2366 

41p3  (a360) 

0.6824 

0.0011 

0.0123 

0.2783 

0.009 

0.0168 

142p4  (a361) 

0.5557 

0.0025 

0.0164 

0.3873 

0.0050 

0.0329 

Table  11:  Summary  of  elemental  scatter  er  distributions  for  Cameron  decomposition  for  Standard 

Quad-Pol  images 


Ship  name 

trihedral 

diplane 

dipole 

cylinder 

narrow 

diplane 

quarter 

wave 

17  Oct.  2005, 122p2  (a342) 

Gulf  Service 

0.2509 

0.0215 

0.0287 

0.2616 

0.0322 

0.4051 

HMSC  Toronto 

0.0544 

0.0544 

0.2646 

0.2840 

0.1790 

0.1634 

20  Oct.  2005,  Dominion  Victory 

141pl  (a358) 

0.0781 

0.0052 

0.2604 

0.3073 

0.1041 

0.2448 

141p3  (a360) 

0.7545 

0.0120 

0.0239 

0.1976 

0 

0.0120 

142p4  (a361) 

0.5849 

0.0018 

0.0055 

0.3894 

0.0018 

0.0166 

6.3.3  SSCM 

The  SSCM  results  are  presented  in  term  of  the  Poincare  sphere  angles,  longitude  2\|/  and  latitude 
2%,  in  Figure  15.  A  ship  is  decomposed  into  six  elemental  scatterers,  as  in  the  Cameron  method. 
A  representative  SSCM  result  presented  in  terms  of  six  types  of  symmetric  scatterers  is  shown  in 
Figure  16  and  the  distribution  of  each  component  is  shown  in  the  histogram.  The  types  of 
elemental  scatterers  are  indicated  as  follows:  1:  trihedral,  2:  diplane  (i.e.  dihedral),  3:  dipole,  4: 
cylinder,  5:  narrow  diplane  and  6:  quarter- wave  device. 

The  results  for  the  cases  considered  are  presented  in  Table  12  and  Table  13. 
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Figure  15:  SSCM  decomposition  images  of  Dominion  Victory  (I41pl,  20  October  2005).  Top:  Fine 
Quad-Pol  Mode;  bottom:  Standard  Quad-Pol  Mode.  Left:  Latitude  coordinate;  right: 

Longitude  coordinate. 

Azimuth  pixels  are  along  the  x-axis  and  range  pixels  are  along  the  y-axis. 
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Figure  16:  SSCM  classification  image  (left)  and  histogram  (right)  of  Dominion  Victory  (I41pl,  20 
October  2005).  Top:  Fine  Quad-Pol  Mode;  bottom:  Standard  Quad-Pol  Mode. 

In  the  histograms  (right):  the  scatter  er  type  is  along  the  x-axis  and  the  percentage  of  the 

distribution  is  along  the  y-axis. 


In  the  SSCM  analysis,  the  assessment  method  of  scattering  coherence  of  the  target  is  different 
from  that  of  the  Cameron  method.  The  scattering  coherence  must  be  tested  before  classification. 
For  the  cases  studied,  the  target  was  tested  as  the  coherent  target  in  only  three  images;  the 
classification  results  are  not  very  promising  since  many  elemental  scatterers  can  not  be  classified. 
This  may  due  to  the  coarse  resolution.  Therefore,  additional  examples  from  a  134  m  long  ship, 
HMCS  Toronto ,  are  given  in  Figure  17  and  Figure  18.  The  results  are  similar  to  those  from  the 
Cameron  method  (Annex  D).  This  method  may  not  be  sufficient  for  small  ship  classification, 
however,  it  can  be  applied  to  large  ships.  The  SSCM  first  finds  the  coherent  scatterers  and  then 
classifies  them.  This  results  in  small  number  of  classified  pixels  of  small  ships  (See.  Figure  15 
and  16).  The  amount  of  information  generated  by  SSCM  in  these  cases  is  not  sufficient  for  ship 
classification.  On  the  other  hand,  for  large  ships,  the  SSCM  produces  an  acceptable  number  of 
classified  pixels  (See  Figure  17  and  18). 
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c)  d) 

Figure  1 7:  SSCM  results  of  HMCS  Toronto  (I22p2,  1 7  October  2005)  of  Fine  Quad-Pol  Mode,  a) 
latitude;  b)  longitude;  c)  classification  and  d)  classification  histogram. 

In  a)  and  b):  azimuth  pixels  are  along  the  x-axis  and  range  pixels  are  along  the  y-axis.  In  d):  the 
scatterer  type  is  along  the  x-axis  and  the  percentage  of  the  distribution  is  along  the  y-axis. 
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Figure  18:  SSCM  results  ofHMCS  Toronto  (I22p2,  17  October  2005)  of  Standard  Quad-Pol  Mode. 

a)  latitude;  b)  longitude;  c)  classification  and  d)  classification  histogram. 

In  a)  and  b):  azimuth  pixels  are  along  the  x-axis  and  range  pixels  are  along  they -axis.  In  d):  the 
scatterer  type  is  along  the  x-axis  and  the  percentage  of  the  distribution  is  along  the  y-axis. 


Table  12:  Summary  of  element  scatterer  distributions  of  SSCM  decomposition  for  Fine  Quad-Pol 

Mode. 


Ship  name 

trihedral 

diplane 

dipole 

cylinder 

narrow 

diplane 

quarter 

wave 

17  Oct.  2005, 122p2  (a342) 

Gulf  Service 

0 

0 

0 

0 

1 

0 

HMSC  Toronto 

0.0598 

0 

0.03424 

0.2120 

0.1848 

0.2011 

20  Oct.  2005,  Dominion  Victory 

141pl  (a358) 

0 

0 

0.6667 

0.2222 

0 

0.1111 
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Table  13:  Summary  of  element  scatter  er  distributions  of  SSCM  decomposition  for  Standard 

Quad-Pol  Mode. 


Ship  name 

trihedral 

diplane 

dipole 

cylinder 

narrow 

diplane 

quarter 

wave 

17  Oct.  2005, 122p2  (a342)  | 

Gulf  Service 

0 

0 

0 

0 

1 

0 

HMSC  Toronto 

0.0213 

0.0106 

0.2872 

0.2128 

0.1809 

0.2872 

20  Oct.  2005,  Dominion  Victory 

141pl  (a358) 

0 

0 

0.75 

0.25 

0 

0 
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7  Summary  and  conclusions 


The  simulated  RADARS AT-2  polarimetric  products  from  the  MARSIE  trial  data  have  been  used 
to  study  the  polarimetric  signature  of  ships.  The  RADARSAT-2  simulator  (RSAT2SIMU) 
simulates  RADARSAT-2  products  using  EC  CV-580  single-look  complex  data  by  increasing  the 
noise  floor  and  degrading  the  spatial  resolution.  However,  the  simulator  can  not  simulate  the 
geometry  of  RADARSAT-2,  such  as  the  altitude,  incidence  angle,  and  swath  coverage.  The 
simulated  data  quality  depends  on  the  input  CV-580  data  quality. 

The  results  demonstrate  that  polarimetric  ship  detection  algorithms  [4]  and  the  classification 
methods,  such  as  target  decomposition  ([5],  [6],  [7]),  developed  using  airborne  PolSAR  data  will 
potentially  be  applicable  to  RADARSAT-2  products. 

A  total  of  5  simulated  images  of  various  ships  have  been  studied  to  date;  it  is  clear  that 
polarimetric  SAR  can  be  used  to  improve  ship  detection  and  to  provide  the  surveillance  operator 
with  some  additional  classification  information.  Polarimetric  SAR  provides  both  phase  and 
amplitude  information,  which  can  reduce  the  false  alarm  rate  and  permit  the  detection  of  smaller 
ships  than  possible  with  single-  or  dual-channel  SAR  systems.  For  example,  the  ship,  Dominion 
Victory ,  having  a  length  of  25  m,  was  not  detected  using  either  single  channel  in  most  cases,  or  by 
a  dual-channel  system  in  some  cases. 

For  the  ship  detection,  the  results  clearly  show  the  advantage  of  polarimetry.  Small  ships  often 
cannot  be  detected  using  single  channel  systems.  The  dual  co-polarization  system  with  phase 
provides  better  detection  performance  than  single  channel  HH  polarization.  However,  dual  co¬ 
polarization  system  should  provide  wider  coverage  than  polarimetric  systems. 

Three  polarimetric  target  decomposition  methods  were  applied  to  characterize  the  targets  of 
interest  in  terms  of  their  elemental  scatterers.  These  methods  included  Pauli,  Cameron,  and 
SSCM.  Each  method  provided  different  information  about  the  target.  The  Cameron  and  SSCM 
target  decomposition  methods  provide  unique  information  on  the  distribution  and  types  of 
scatterers  that  constitute  the  target  of  interest,  providing  an  approach  for  target  classification. 
These  decomposition  methods  work  best  for  well-focused,  large  ships.  The  Pauli  decomposition 
method  provides  a  simple  but  efficient  way  to  distinguish  the  target  of  interest  (predominantly 
double  bounce  scattering)  from  the  surrounding  ocean  (predominantly  single  bounce  scattering). 

Decomposition  results  from  the  considered  methods  should  be  combined  with  other  target 
metreces  such  as  signature  length,  signature  cross  section,  and  derived  tar  motion  to  aid  in  target 
classification. 

Based  upon  our  analysis,  the  following  issues  are  recommended  for  further  study: 

•  The  relationship  between  target  features  and  the  elemental  scatterers  derived  from  the 
Cameron  and  the  SSCM  polarimetric  decomposition  methods; 

•  The  relationship  between  target  features  and  the  polarimetric  signature;  and 
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•  The  feasibility  of  applying  polarimetric  target  analysis  methods  to  automatic  target 
recognition  by  using  scattering  elements  as  target  features  in  target  classification  software 
such  as  the  ATR  Workbench  [15]. 
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Annex  A  Simulated  RADARSAT-2  images 


The  simulated  RADARSAT-2  images  of  Dominion  Victory  from  four  flight  passes  on  20  October 
2005  and  one  flight  pass  which  contains  three  known  vessels  on  17  October  2005  are  presented  in 
this  Annex.  The  Dominion  Victory  from  141p2  and  Edward  CornWallis  from  122p2  are  not 
visible  in  the  simulated  images. 

In  the  figures  in  this  Annex,  azimuth  pixels  are  along  the  x-axis  and  range  pixels  are  along  the  y- 
axis. 


100  200  300  400  500  600 


a) 


b)  c) 

Figure  19:  Polarimetric  image  of  Dominion  Victory  (I41p2,  20  October  2005).  a)  CV-580  image; 

b)  Fine  Quad-Pol  Mode;  c)  Standard  Quad-Pol  Mode. 
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b)  c) 

Figure  20:  Polarimetric  image  of  Dominion  Victory  (I41p3,  20  October  2005).  a)  CV-580  image; 

b)  Fine  Quad-Pol  Mode;  c)  Standard  Quad-Pol  Mode. 
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Figure  21:  Polarimetric  image  of  Dominion  Victory  (I42p4,  20  October  2005).  a)  CV-580  image; 

b)  Fine  Quad-Pol  Mode;  c)  Standard  Quad-Pol  Mode. 
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b) 


c) 


Figure  22:  Polarimetric  image  of  Gulf  Service  (I22p2,  1 7  October  2005).  a)  CV-580  image;  b) 

Fine  Quad-Pol  Mode;  c)  Standard  Quad-Pol  Mode. 


38 


DRDC  Ottawa  TM  2007-189 


b) 


c) 


Figure  23:  Polarimetric  image  of  Edward  Cornwallis  (I22p2,  1 7  October  2005).  a)  CV-580 
image;  b)  Fine  Quad-Pol  Mode;  c)  Standard  Quad-Pol  Mode. 


DRDC  Ottawa  TM  2007-189 


39 


100  200  300  400  500  600  700  800  900  1000  1100 


b)  c) 

Figure  24:  Polarimetric  image  ofHMCS  Toronto  (I22p2,  17  October  2005).  a)  CV-580  image;  b) 

Fine  Quad-Pol  Mode;  c)  Standard  Quad-Pol  Mode. 
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Probability  of  missed  detection 


Annex  B  ROC 


The  ROC  results  of  cases  studied  are  presented  in  this  Annex. 


Figure  25:  Detection  performance  for  Dominion  Victory  (I41p2,  20  October  2005).  a)  Fine  Quad-Pol 

Mode;  b).  Standard  Quad-Pol  Mode. 


a  b 

Figure  26:  Detection  performance  for  Dominion  Victory  (I41p3,  20  October  2005).  a).  Fine 

Quad-Pol  Mode;  b).  Standard  Quad-Pol  Mode. 
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Probability  of  missed  detection 


Figure  27:  Detection  performance  for  Dominion  Victory  (I42p4,  20  October  2005).  a)  Fine  Quad- 

Pol  Mode;  b)  Standard  Quad-Pol  Mode. 
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Annex  C  Pauli  Decomposition  Results 


The  Pauli  decomposition  results  for  the  cases  studied  are  presented  in  this  Annex. 


Figure  28:  Pauli  decomposition  image  of  Dominion  Victory  (I41p3,  20  October  2005).  a)  Fine 

Quad-Pol  Mode;  b)  Standard  Quad-Pol  Mode. 
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Figure  29:  Pauli  decomposition  image  of  Dominion  Victory  (I42p2,  20  October  2005).  a)  Fine 

Quad-Pol  Mode;  b)  Standard  Quad-Pol  Mode. 
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Figure  30:  Pauli  decomposition  image  of  Gulf  Service  (I42p2,  17  October  2005).  a)  Fine  Quad- 

Pol  Mode;  b)  Standard  Quad-Pol  Mode. 


120  140  160  180  200  220  240 

Azimuth  Pixels 


120  140  160  180  200  220  240 

Azimuth  Pixels 


a)  b) 

Figure  31:  Pauli  decomposition  image  ofHMCS  Toronto  (I22p2,  17  October  2005).  a)  Fine 

Quad-Pol  Mode;  b)  standard  quad-plo  Mode. 
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Annex  D  Camerom  decomposition  results 


The  Cameron  decomposition  results  of  cases  studied  are  presented  in  this  Annex. 

In  the  histogram  figures  in  this  Annex:  the  scatterer  type  is  along  the  x-axis  and  the  percentage  of 
the  distribution  is  along  the  y-axis. 
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Figure  32:  Cameron  decomposition  image  and  histogram  of  Dominion  Victory(  I41p3,  20 
October  2005).  a)  Fine  Quad-Pol  Mode;  b)  Standard  Quad-Pol  Mode. 
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Figure  33:  Cameron  decomposition  image  and  histogram  of  Dominion  Victory  (I42p4,  20 
October  2005).  a)  Fine  Quad-Pol  Mode;  b)  Standard  Quad-Pol  Mode. 
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Figure  34:  Cameron  decomposition  image  and  histogram  of  Gulf  Servic  (I22p2,  1 7  October 
2005).  a)  Fine  Quad-Pol  Mode;  b)  standard  quad-plo  Mode. 
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Range  Pixels  Range  Pixels 
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Figure  35:  Cameron  decomposition  image  and  histogram  of  HMCS  Toronto  (I22p2,  1 7  October 
2005).  a)  Fine  Quad-Pol  Mode;  b) standard  quad-plo  Mode. 
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Annex  E  SSCM  decomposition  results 


The  results  of  the  SSCM  decomposition  are  presented  in  this  Annex.  There  weren’t  any  coherent 
targets  detected  in  141p3  or  142p4. 

In  the  figures  in  this  Annex,  a)  and  b):  azimuth  pixels  are  along  the  x-axis  and  range  pixels  are 
along  the  y-axis;  d):  the  scatterer  type  is  along  the  x-axis  and  the  percentage  of  the  distribution  is 
along  the  y-axis. 
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Figure  36:  SSCM  decomposition  images  Gulf  Service  for  Fine  Quad-Pol  Mode  (I22p2,  17 
October  2005).  a)  latitude;  b)  longitude;  c)  classification;  d)  histogram. 
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Figure  3  7:  SSCM  decomposition  images  of  Gulf  Service  for  Standare  Quad-Pol  Mode  (I22p2,  1 7 
October  2005).  a)  latitude;  b)  longitude;  c)  classification;  d)  histogram. 
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Figure  38:  SSCM  decomposition  images  of  Toronto  for  Fine  Quad-Pol  mode  (I22p2,  17  October 
2005).  a)  latitude;  b)  longitude;  c)  classification;  d)  histogram. 
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Figure  39:  SSCM  decomposition  images  of  Toronto  for  Standard  Quad-Pol  mode  (I22p2,  17 
October  2005).  a)  latitude;  b)  longitude;  c)  classification;  d)  histogram. 
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Annex  F  Ship  photos 


The  photographs  of  the  ships  studied  are  shown  in  this  Annex.  These  photos  were  taken  during 
the  MARSIE  trial  by  Janice  Lang  (DRDC  Ottawa). 


a)  b) 


c) 

Figure  40:  Ship  photographs:  a)  Dominion  Victory;  b)  HMCS  Toronto;  c)  Gulf  Service. 
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XML 


Canada  Center  for  Remote  Sensing 
Canadian  Space  Agency 
Chip-Based  Adaptive  SAR  Processor 
Configurable  Airborne  SAR  Processor 
Defence  R&D  Canada 
Environment  Canada 
Geographic  Tagged  Image  File  Format 
Horizontal-Horizontal 
Horizontal- V  ertical 
Look-Up  Table 

Maritime  Sensor  Integration  Experiment 
Noise  Equivalent  Signal  Zero 
Probability  of  False  Alarm 
Probability  of  Missed  Detection 
Peak-to-Clatter  Ratio 
Polarimetric  SAR  Interferometry 
Polarimetric  Synthetic  Aperture  Radar 
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Single  Look  Complex 
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